Chirality dependent frequency shift of radial breathing mode in metallic carbon 
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A phonon frequency shift of the radial breathing mode for metallic single wall carbon nanotubes is 
predicted as a function of Fermi energy. Armchair nanotubes do not show any frequency shift while 
zigzag nanotubes exhibit phonon softening, but this softening is not associated with the broadening. 
This chirality dependence originates from a curvature-induced energy gap and a special electron- 
phonon coupling mechanism for radial breathing modes. Because of the particle-hole symmetry, only 
the off-site deformation potential contributes to the frequency shift. On the other hand, the on-site 
potential contributes to the Raman intensity, and the radial breathing mode intensity is stronger 
than that of the G band. The relationship between the chirality dependence of the frequency shift 
of the radial breathing mode and the F point optical phonon frequency shift is discussed. 
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I. INTRODUCTION 

Raman spectroscopy has been widely used for 
the characterization of carbon nanotubei^^^ and 
graphene 7 ^! 9 - samples, since Raman spectroscopy is a 
non-destructive and non-contact measurement that can 
be carried out under ambient air pressure and at room 
temperature. The G band which gives a strong Raman 
intensity at around 1600cm -1 (~0.2eV) in graphene con- 
sists of the r point longitudinal optical (LO) and trans- 
verse optical (TO) phonon modest For single wall car- 
bon nanotubes (SWNTs), these two phonon modes with 
A\ g symmetry split by the curvature effect into two Ra- 
man features which we call the G + and G~ bands^i It 
has been observed for metallic SWNTs that the peak po- 
sitions of the G + and G~ features shift differently from 
each other as a function of the Fermi energy, j^pi 12 ! 13 ! 14 ! 15 
It is known that virtual electron-hole pair creation by the 
electron-phonon (el-ph) interaction gives a self-energy 
correction to the phonon frequency, which is relevant to 
the phonon frequency shift. The behavior of the G + 
and G~ bands indicates that the el-ph interactions for 
the r point LO/TO modes appear to be different from 
each other in metallic SWNTs. The frequency shift for 
the LO/TO mode for graphene and nanotubes has been 
discussed by many authors, 16 . 17 i 18 i l9 i 20 . 21 i 22 . 23 i 24 i 25 In a 
previous paper, we showed that the phonon softening for 
the LO and TO phonon is chirality dependent in which 
the curvature effect is important^ For example, only 
armchair SWNTs do not exhibit any frequency shift of 
the TO mode due to the absence of a curvature-induced 
energy gap in armchair SWNTs. 

A similar effect may be observed for the radial breath- 
ing mode (RBM) which appears in the SWNT Raman 
spectra at around 240cm~ 1 (~30meV) for a nanotubc 
diameter of around d t ~ lnm. The RBM is often used to 
assign the diameter and chirality of a nanotube^ i 27 ' 28 For 
metallic SWNTs except for armchair SWNTs, the curva- 



ture of a cylindrical graphene layer induces a small energy 
gap. 29 ' 30 ! 31 Since the curvature-induced energy gap has a 
similar energy to that for the RBM, especially for SWNTs 
with diameter less than 2nm, the curvature-induced en- 
ergy gap may affect the phonon softening of the RBM, 
which is the motivation of the present paper. Thus, it is 
important to estimate this frequency shift of the RBM 
for metallic SWNTs. In this paper, we examine the de- 
pendence of the RBM frequency shift on E-p and chirality 
for metallic SWNTs. The relationship between the RBM 
and the LO/TO modes is examined for the el-ph matrix 
element of the electron-hole pair creation. 

Since the frequency shift of a phonon mode is mainly 
due to a low energy electron-hole pair creation, a theory 
for electrons near the Fermi energy is useful for calculat- 
ing the matrix element for an electron-hole pair creation. 
The el-ph interaction has previously been examined in 
terms of an effective-mass model by several authorsi 19 i 32 
It is known that the el-ph interaction consists of on-site 
and off-site terms. The on-site term represents a process 
whereby a 7r-electron changes its energy due to a defor- 
mation potential but stays at the same position. The 
off-site term on the other hand represents a scattering 
process whereby a 7r-electron moves into a nearest carbon 
site due to a deformation potential^ As for the T point 
LO/TO modes, the el-ph interaction is given only by the 
off-site term. In a previous paper, we showed that the off- 
site term of the el-ph interaction becomes sensitive to the 
position of the cutting line or to the curvature effect^ In 
this paper, we will show that an interesting property of 
the RBM is that the el-ph interaction is given not only by 
the off-site term but also by the on-site term. Although 
the on-site term does not contribute to the frequency shift 
of the RBM because of the particle-hole symmetry about 
the Dirac point, the on-site term of the el-ph interaction 
enhances the Raman intensity of the RBM. 

This paper is organized as follows. In Sec.|TT]we show 
that the curvature-induced energy gap determines the 
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basic features for the frequency shift of the RBM. In 
Sec. IIIIl we show the el-ph interaction for the RBM using 
an effective-mass model for a graphene sheet with a lat- 
tice deformation. In Sec. IIV1 we calculate the frequency 
shift of the RBM as a function of the Fermi energy. We 
compare our results with experimental data in Sec. [V] 
A discussion and summary of these effects are given in 
Sec. IVI1 The effective-mass Hamiltonian used in Sec. IIIIl 
and Sec. IIVI is derived in Appendix [Al and Appendix [Bl 
respectively. 



II. RADIAL BREATHING MODE AND THE 
CURVATURE-INDUCED ENERGY GAP 

Here we show that the basic features of the frequency 
shift of the RBM are determined by the curvature- 
induced energy gap. A renormalized RBM energy be- 
comes Tito = fiLj( ' + Huj^ where u>^ is the unperturbed 
RBM frequency and a/ 2 -* is the quantum correction to the 
RBM frequency due to the el-ph interaction. We assume 
that is a monotonic function of the tube diameter 
(dt [ran] ) and is modeled as linear in inverse diameter with 
an offset as 



ci 
d t 



c 2 , 



(1) 



where c\ = 223.5[cm _1 ] and c 2 = 12.5[cm _1 ] are exper- 
imentally derived parameters as obtained by Bachilo et 
aL ) 34 ' 35 while Huj 1 - 2 ^ is calculated by second-order pertur- 
bation theory 2 ^ as 



|<eh(k)|W, 



cp 



,(0)\ 



7kj(°) - [E c (k) - E h (k)]+iT 



{f[E h (k)-E F ]-f[E e (k)-E F }}. 



(2) 



In Eq. ([2]), the factor 2 comes from spin degeneracy, 
E c (k) (i?h(k)) is the energy of an electron (hole) with 
wave vector k, (eh(k)|7i op |u; ( ' ' ) ) is the el-ph matrix ele- 
ment that the RBM changes into an electron-hole pair 
with wave vector k (see Fig. [Ha)), and f[x] is the Fermi 
function. We obtain T (the half of the decay width) by 
calculating T = — Im(?ia/ 2 ') self-consistently in Eq. ((2|). 

We consider the real part and imaginary part of the 
denominator of Eq. ©, h(E) = l/(7kJ™ - E + iT), 
as a function of the energy of an intermediate electron- 
hole pair state, E = E e (k) ~ E h {k) (see Fig. Hfb)). 
Ke(h(E)) has a positive (negative) value when E < huj^ 
(E > huj^), and the lower (higher) energy electron-hole 
pair makes a positive (negative) contribution to huj( 2 \ 
Therefore, the sign of the contribution to Re(?iu/ 2 )), 
i.e., frequency hardening or softening, depends on its 
electron-hole virtual state energy, E. The curvature- 
induced energy gap, E gap , affects the RBM frequency 
shift since an electron-hole pair creation is possible only 
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FIG. 1: (Color online) (a) An intermediate electron-hole pair 
consists of an electron in the conduction band and a hole in 
the valence band. The RBM is denoted by a zigzag line and 
an electron-hole pair is represented by a loop. A cutting line 
for a metallic SWNT is shown as the solid curve projected on 
the surface of the Dirac cone. The electron-hole pair creation 
is possible only when E > E ga _ p . (b) The energy correction 
to the phonon energy by an intermediate electron-hole pair, 
especially the sign of Ke(h(E)) (black curve), that is, fre- 
quency hardening or softening, depends on the energy of the 
intermediate state E. The contribution to hui^ of a low en- 
ergy electron-hole pair satisfying < E < iJ gap is forbidden. 
lm(h(E)) (red curve) is nonzero only when E is very close to 
Huj( \ which shows that the RBM can resonantly decay into 
an electron-hole pair with the same energy. The figure is the 
case of huj m = 30meV, £ gap = 20meV. The Fermi distribu- 
tion function /[_Eh(k) — E-p] — f[E c (k) — Ef] is plotted for 
T — 10K and T = 300K as a function of E in the case of 
\E F \ = 20meV. 



When E gap > huj^°\ not only are all the positive con- 
tributions to the RBM frequency suppressed, but some 
negative contributions are also suppressed. Further, 
Im(h(E)) is nonzero only when E is very close to Tiuj^ \ 
which shows that the RBM phonon can resonantly de- 
cay into an electron-hole pair with the same energy. 
This means that T ~ when E gap > Tiuj^ because no 
electron-hole pair excitation is allowed near E = huj(°\ 



bution to frequency hardening in Eq. (J2]) is suppressed. Thus, it is important to compare the values of E gap and 
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huW for each (n,m) SWNT. 
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d t (nm) 



fkjW w hen d t < 2nm (see Fig. E^b)). Thus, no de- 
cay process (or spectral broadening) contributing to T 
is available for a zigzag SWNT with d t < 2nm. On the 
other hand, a decay process contributing to T is avail- 
able regardless of the chirality of a metallic SWNT with 
dt > 2nm. For example, in Fig.[2jb), we see that a (12, 0) 
tube does not exhibit a finite value of T even though we 
get a finite value for the el-ph matrix element, but (9, 6) 
and (10, 7) may exhibit a broadening. 

In addition, because of the Fermi distribution function 
f[x] in Eq. an electron-hole pair satisfying E < 2\Ep\ 
can not contribute to the energy shift at zero tempera- 
ture due to the Pauli principle. At zero temperature, 
when |E F | ~ hu^ /2, then hui^ takes a minimum value 
since all positive contributions to huj^ are suppressed in 
Eq. ([2]). 26 Ep and E gap play a very similar role at zero 
temperature, but their difference becomes clear at a fi- 
nite temperature. For example, when E gap < ftio^\ then 
F can be nonzero even when 2|Ep| > hu)^ '. On the other 
hand, when -Egap > huj^°\ then Tius^ does not have an 
imaginary part (T = 0) even at room temperature, re- 
gardless of the value of E-p. This difference between Ep 
and i?g a p is understood in Eq. @ as (eh(k)|W ep |a/°)) = 
for E < Egap and as that f[E h (k) - Ep] - f[E c {k) - Ep] 
is not zero for E < 2\Ep\ (at a finite temperature) as 
shown in Fig. [ljb). Thus, E gap determines whether a 
SWNT can exhibit a broadening. 



III. ELECTRON-PHONON INTERACTION BY 
EFFECTIVE-MASS THEORY 



FIG. 2: (Color online) (a) The dependence of the curvature- 
induced energy gap, -E gap , on the chiral angle 6 and tube di- 
ameter d t . The surface is a plot of Eq. © which reproduces 
well the calculated results, (b) The d t dependence of -E ga p, 
which is given as a one-dimensional projection of (a) onto the 
dt axis. We plot the energy of the RBM, 7kj (0) of Eq. (JTJ , as 
a solid blue curve for comparison. The points on the dashed, 
dot-dashed, and dotted curves satisfy n — m = 3, 6, 9, respec- 
tively. 

In Fig. [2{a) we plot E gap for each (n, m) for metallic 
SWNTs as a function of the chiral angle 0(°) and tube di- 
ameter dt (nm). We performed the energy band structure 
calculation in an extended tight-binding framework 3 ^ to 
obtain E gap . Figure [2]^a) shows that, for a fixed dt, a 
zigzag SWNT (0 = 0°) has the largest value of E gap 
and an armchair SWNT (9 = 30°) has no energy gap. 
The (n, m) values associated with a curve are given by 
n — m = 3q (q = 0, 1, 2, • • • ). The calculated results are 
well reproduced by 



Egap = cos 39, 



with c = 60(meV-nm 2 ) 
zigzag SWNTs (0 = 0°), 



M- Using Eqs. 
we see that E £ 



(3) 

© and © for 
ap is larger than 



Next we show the el-ph interaction in the Hamiltonian 
by effective-mass theory, which is used for calculating the 
matrix element of the el-ph interaction. We will derive 
Eqs. (g]), (JHJ) and © from the nearest-neighbor tight- 
binding Hamiltonian in Appendix [Al 38 



A. Unperturbed Hamiltonian 

The unperturbed Hamiltonian in the effective-mass 
model for 7r-electrons near the K point of a graphcnc 
sheet is given by 



H Q = VpCT ■ p, 



(4) 



where vp is the Fermi velocity, p = — ih\7 is the momen- 
tum operator, and cr = (cr x , cr y ) is the Pauli matrix. Ti.^ 
is a 2 x 2 matrix which operates on the two component 
wavefunction: 



^ K (r) 



^(r) 
V£(r) 



(5) 



where V'aI 1 ') 



and ipftir) are the wavefunctions of ir- 
electrons for the sublattices A and B, respectively, around 
the K point. The energy eigenvalue of Eq. Q is given 
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by ±uf|p| and the energy dispersion relation shows a lin- 
ear dependence at the Fermi point ^ which is known as 
the Dirac cone. The eigenstates for E = +uf|p| and 
E = — fF|p| are given by 



(a) 



V< c K k (r) = 
< k (r) = 



e -i0(k)/2 

e + 4 e(k)/2 

e -iS(k)/2 

_ e +ie(k)/2 



(6) 



which are a conduction state and a valence state with 
wavevector k, respectively. In Eq. ©, S denotes the 
surface area of graphene, the wavevector is measured 
from the K point, and O(k) is defined by an angle of 
k = (k x ,ky) measured from the k x axis as (k x ,k y ) = 
|k|(cos0(k),sinO(k)). Here the k y axis is defined by 
the direction of a zigzag nanotube axis (see the coor- 
dinate system in Fig. [3)Ja)). The energy eigenstate for 
the valence band, V^k( r ) ^ s gi ven by a z ip^ k (r). This 
results from the particle-hole symmetry of the Hamilto- 
nian: a z Ti^a z = —Tig. The dynamics of 7r-electrons 
near the K' point relates to the electrons near the K 
point by time- reversal symmetry, ip K — > (ip K )*. Because 
lattice vibrations do not break time-reversal symmetry, 
we only consider the electrons near the K point in this 
paper. 



B. Perturbation 

Lattice deformation modifies the nearest-neighbor hop- 
ping integral locally as —70 — > —70 + Sj§ (r,) (a = 1,2,3) 
(see Fig. [3](a)). The corresponding perturbation of the 
lattice deformation is given by 



jgA o=l,2,3 



where cf" is the annihilation operator of a ir electron of an 
A- atom at position r^, and (cf +a Y is a creation operator 
of position r^ +a (= + R a ) of a B-atom where R a (a = 
1,2,3) are vectors pointing to the three nearest-neighbor 
B sites from an A site (see Fig. [3ta)). 

The perturbation of Eq. ([7]) gives rise to scattering 
within a region near the K point (intravalley scatter- 
ing) whose interaction is given by a deformation-induced 
gauge field A(r) = (A x (r), A y (r)) in Eq. gj as 



vpa ■ [p + A(r)] . 
A(r) is defined from Sjq(t) (a = 1,2,3) a; 



(8) 




(b) s,(x)n(xl , s ( a .) 




*s~(x + dx)n(x + dx) 
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FIG. 3: (a) A hexagonal unit cell of graphene consists of 
A (closed circle) and B (open circle) atoms, ai and a.2 are 
lattice vectors. R a (a = 1,2,3) are vectors pointing to the 
nearest-neighbor B sites from an A site. For the coordinate 
system of (x, y), the R a (a = 1, 2, 3) are written, respectively, 
as Ri = a cc e y , R 2 = — (v / 3/2)a cc e ;c — (l/2)a cc e y , and R 3 
(V3/2)a cc e x 



(l/2)a cc e H where a cc is the carbon-carbon bond 



length and e x (e a ) is the dimensionless unit vector for the x- 
axis (y-axis). Local modulations of the hopping integral are 
defined by <5 7 o(ri) (a = 1,2,3) where ri is the position of 
an A-atom. (b) The displacement vector for the RBM, s(x), 
is decomposed in terms of the normal s z (x) and tangential 
s x (x) components. The derivative of the normal unit vector 
n(x) with respect to x gives a component along e x (x), which 
modifies the net displacement along the x direction. 



A(r) • R2 = 0. Generally, the direction of A(r) is point- 
ing perpendicular to the bond whose hopping integral is 
changed from 70. 

When the displacement vector of a carbon atom at 
r is s(r) = (s x (r), s y (r), s z (r j), the perturbation to the 
nearest-neighbor hopping integral is given by 



5 7o °(r) 



ffoff 

£(lc.c 



R a -{s(r + R a ) -s(r)}, 



(10) 



,38,39 



where g s is the off-site coupling constant and t — 3a cc /2. 
By expanding s(r + R Q ) in a Taylor series around s(r) as 
s(r + R a ) = s(r) + (R a ■ V)s(r) + • • • , we approximate 
Eq. (£TDJ) as 



v F A x (r) = SjKt) - - [<5 7o 2 (r) + S^r)} , 
v F A y (r) = ^[8 7 *(r)-5 7 3 (r)]. 



(9) 



When 5^1 = 5^ = 0, then A(r) = (A x (r),0) and 
A(r) • Ri = 0. Similarly, when = <57q = 0, we have 



<5 7o Q (r) * |^R Q ■ {(R a • V)s(r)} . (11) 



Putting Ri = a cc e y , R2 = -(V3/2)a cc e x - (l/2)a cc e y , 
and R 3 = (yZ/2)a cc e x — (l/2)a cc e y , into the right- 
hand side of Eq. (|lip . we obtain the corresponding 
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deformation-induced gauge field of Eq. ([9]) as 



vfA x (t) 
v F A y (r) 



goff 

2 

5off 



d^(r) _,_ 9s y (r) 
dx 
ds x (r) 



dy 



dy 
ds y (r) ~ 

dx 



(12) 



Further the displacements of carbon atoms give an on- 
site deformation potential 



'Ho 



ds x (r) dsyjr) 
dx dy 



(13) 



Here Oo is the 2x2 unit matrix and ds x (r)/dx + 
ds y (r) /dy represents the change of the area of a graphene 
sheet. 32 In Eqs. (12) and (13), according to the density 
functional calculation by Porezag et aZ.j^S we adopt the 
off-site coupling constant g Q g = 6.4eV and the on-site 
coupling constant g on = 17.0eVi 33 i 41 

Since Eqs. (fl~2|) and (fl"3|) are proportional to the deriva- 
tives of s x (r) and s y (r), that is, they are proportional 
to q, the el-ph matrix element for the in-plane longitu- 
dinal/transverse acoustic (LA/TA) phonon modes van- 
ishes at the r point where q is the phonon wave vector. 
Namely, A(r) = and H on = in the limit of q = 0. 
Among the TA phonon modes, there is an out-of-plane 
TA (oTA) phonon mode. The oTA mode shifts carbon 
atoms on the flat 2D graphene sheet into the z-direction 
(see Fig. |3][a) and (b)). The oTA mode of graphene cor- 
responds to the RBM of a nanotube even though the 
RBM is not an acoustic phonon mode. 10 In the follow- 
ing, we will show that the el-ph interaction for the RBM 
is enhanced due to the curvature of the nanotube as com- 
pared with the oTA mode of graphene since the RBM is 
a bond-stretching mode because of the cylindrical struc- 
ture of SWNTs. 

The displacements of the RBM modify the radius of 
a nanotube as r — > r + s z (r) (see Fig. [3]Jb)). A change 
of the radius gives rise to two effects on the electronic 
state. One effect is a shift of the wavevector around the 
tube axis. The distance between two wavevectors around 
the tube axis depends on the inverse of the radius due 
to the periodic boundary condition and a change of the 
radius results in a shift of the wavevector. The other 
effect is that the RBM can change the area on the surface 
of the nanotube even at the T point. This results in 
an enhancement of the on-site interaction. These two 
effects are relevant to the fact that the normal vector 
on the surface of a nanotube is pointing in a different 
direction depending on the position. To show this, we 
take a (zigzag) nanotube as shown in Fig. \3[b). Let us 
denote the displacement vectors of two carbon atoms at 
x and x + dx as s(x) and s(x + dx), then an effective 
length for the displacement along the x axis between the 
nearest two atoms is given by 



D x = e x (x + dx) ■ [s(x + dx) — s(x)] 



(14) 



By decomposing s(x) in terms of a normal and a tan- 
gential unit vector as s(x) = s z (x)n(x) + s x (x)e x (x) (see 



Fig-EJb)), we see that Eq. (jT4j) becomes 

D x = s x (x + dx) + s z (x + dx)e x (x + dx) ■ n(x + dx) 
— s x (x)e x (x + dx) ■ e x (x) — s z (x)e x (x + dx) ■ n(x) 
ds x (x) s z (x) 



dx 



dx 



(15) 



where we have used the following equations: 

n(x + dx) = n(x) + —e x (x) H 

r 

e x (x + dx) = e x (x) + —nix) H 

r 



(16) 



Equation (fl~5|) shows that the net displacement along the 
x axis is modified by the curvature of the nanotube as 
d x s x (r) — > d x s x (r) + s z (r) jr. The correction is negligible 
for a graphene sheet (r — > oo). 

The el-ph interaction for the RBM is included by re- 
placing d x s x (r) with d x s x (r) + s z (r)/r in Eqs. (fl"2")) and 
(p~3|) . In Eq. (| 1 2[> . we have an additional deformation- 
induced gauge field, V F A x (r) = -(g oS /2)(s z (r)/r), for 
the RBM mode which gives rise to a shift of the wavevec- 
tor around the tube axis. In Eq. (|13[) . it is shown 
that the RBM produces on-site deformation potential of 
5onO"o(s 2 (r)/r). Finally, we obtain the el-ph interaction 
for the r point (q = 0: s(r) is a constant) RBM, as 



~j 9oS s z s z 

Hep = ~ &x + Son C 



2s 2 



Son 
I 

"2 



<I:H 
2 



(17) 



IV. 



ELECTRON-PHONON MATRIX ELEMENT 
AND FREQUENCY SHIFT 



In this section we calculate the el-ph matrix element 
for an electron-hole pair creation and the corresponding 
frequency shift of the RBM. 

From Eqs. ([6]) and (fl"7]) . the el-ph matrix element for 
an electron-hole pair generation near the K point is given 
by 



(eh(k)|W„ 



,(o)\ 



= / (^ c K k (r))tH cp < k (r)d 2 r 




a 

yon 2 

_9ott a 

2 yon 



.-<->(k) 



"Ik) 



igoff-r sin0(k). 
d t 



(18) 



It is noted that the g on term does not contribute to the 
RBM frequency shift. This is because the particle- hole 
symmetry: Vv,k( r ) = c^c k(r), gives a vanishing matrix 
element, ^l k (r)g on ao^ V:]i (r) = g on V'c,k( r ) cr oc z '0c,k(r) = 
in Eq. (JT8J). We note that sin9(k) of Eq. |TS]) indi- 
cates that low energy electron-hole pairs near the Dirac 
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point (k-states satisfying 8(k) w on the cutting line, 
see Fig. @[a)) are hardly excited. Instead, high energy 
electron-hole pairs (k-states of 0(k) « ±tt/2) do con- 
tribute to the frequency softening. Putting Eq. (jT5J) into 
Eq. ([2]), we calculate the frequency shift as a function 
of E F for a (9,0) zigzag SWNT. In Fig. HTj), we plot 
(black line) and uj for room temperature (red curve) 
and u> for 10 K (blue curve). The frequency difference 
between Ep = 0.6 eV and the Dirac point (-Bp = 0) is 
about 10 cm -1 . 



(a) 




vf\Ic\ = u 



W 0) 




(b) 



(9,0) 



Cutting line 

£ gap rj 0.12(eV) 




-0.6 -0.4 -0.2 0.2 0.4 0.6 

E F (eV) 



FIG. 4: (Color online) (a) The decay width is zero because 
of -Egap > Tuo ° ■ This is shown by the relative positions of 
the cutting line and the equi-energy contour circle satisfying 
vf\^\ = with respect to the Dirac point, (b) The Ef 
dependence of the RBM frequency in the case of the (9, 0) 
zigzag SWNT at room temperature (red curve) and at 10 K 
(blue curve). We plot as the black line for comparison. 

It is useful to compare Eq. I|18p with the amplitude of 
electron- hole pair creation by the T point LO/TO phonon 
modes in order to understand the diameter dependence of 
the RBM frequency shift of a zigzag SWNT. In a previous 
paper we obtained for zigzag SWNTs that 



(eh(k)|H op Ko) 
(eh(k)|H C p|u;To) 



-»5off smB(k), 

^CC 

-«5off cosO(k), 



(19) 



where wlo ( u to) is the amplitude of the LO (TO) 
phonon mode. The O(k) dependence of the matrix el- 
ement of Eq. (fT8|) is the same as the T point LO mode 
of Eq. (TT1]). By comparing Eq. (TS]) with Eq. we 
find that the ratio of the el-ph matrix element squared is 
given by 



R = 



(eh(k)|W 



cp 



,(0)\l 2 



|(eh(k)|W 0p K o )| 2 \d t u LO J ■ (20) 

Since the phonon amplitude is proportional to the 
phonon frequency as s z cx yjTi/Mto^ and ulo oc 
^/S,/A/wlo (M is the mass of the carbon atom), we ob- 
tain from Eq. (J20J) that 



d t 



(21) 



Here, we have used Eq. {1} and wlo = 1600cm -1 to get 
s 1/ u lo ~ w lo/ w ' ' ) — 7.2(d t /lnm). Because the squared 
matrix element appears in the numerator of Eq. (|2|) , the 
RBM frequency shift is proportional to the inverse of 
d t . This is consistent with the fact that the RBM fre- 
quency shift of a (9,0) SWNT is around 10cm -1 because 
the frequency shift of the T point LO mode 2 & reaches 
around 50cm -1 and R ~ 0.2. The frequency difference 
between E-p = 0.6eV and the Dirac point is about 5 cm -1 
for a (18,0) zigzag SWNT (d t ~ 1.4nm), whose d t is 
twice the d t of a (9,0) SWNT. It is expected that the fre- 
quency shift of the RBM in zigzag SWNTs is expressed 
by 50(a cc /d t )cm -1 . 

In Fig. 0Jb) , the decay width is zero because there is 
no electron-hole pair satisfying E < Tav^ ' (see Figs. Hta) 
and[2^b)). In principle, the decay width due to the el-ph 
interaction is absent for zigzag SWNTs when dt ^ 2nm 
because the curvature-induced energy gap is larger than 
the original RBM phonon energy 

: Bgap > ?iw (0) - On the 
other hand, when dt > 2nm, we have checked that the T 
value is less than 1cm -1 due to the small matrix clement 
for a larger diameter zigzag SWNT. 

For a general (n, m) SWNT with a chiral angle 8, the 
el-ph interaction for the RBM becomes 



H cp {9) 



2s, 



ffon 
9off „-i38 



gait r +l38 
2 e 



(22) 



which is derived in Appendix [B] As a result, the ma- 
trix element for an electron-hole pair creation is chirality 
dependent as 

(eh(k)|H cp (#)| = ig oS £ sin(9(k) + 30), (23) 

where ki — |k| cos O(k) is the wavevector in the direction 
around the tube axis and hi = |k|sinO(k) is the one 
along the tube axis. Thus the frequency shift of the RBM 
has a chiral angle dependence. In particular, armchair 
SWNTs (0 = 30°) do not exhibit any frequency shift 
because the matrix element becomes 



(eh(k)|«ep(0)|w 



(0)\ 



igozf cos(8(k)), (24) 
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which is zero for a cutting line of the metallic band: 
6(k) = ±7r/2. This 6(k) dependence of Eq. (H} is the 
same as that of the TO phonon mode of Eq. (|19| and the 
absence of the frequency shift of the RBM in armchair 
SWNTs is similar to that of the T point TO mode in 
armchair SWNTs^ 

In Fig. [H^a) we show the frequency shift in a (9, 6) 
SWNT (9 = 24°). The factor of sin(0(k)+36>) in Eq. QZ® 
indicates that low energy electron-hole pairs satisfying 
O(k) ps 7r/2 — W = 18° on the cutting line contribute 
significantly to the frequency shift. The decay width 
is non-zero because an electron-hole pair can be excited 
near E ~ Tiuj^ (see Figs. EJb) and^b)). 



(a) 




(b) 



-0.15 -0.1 -0.05 0.05 0.1 0.15 

E F (eV) 




up|k| = ^ 0) 




R-2 

Cutting line 



FIG. 5: (Color online) (a) The E F dependence of the RBM 
frequency in the case of the (9, 6) chiral SWNT at room tem- 
perature (red curve) and at 10 K (blue curve). We plot a/ ' 
as the black line for comparison, (b) The decay width (T) is 
non-zero because of -E gap < Tilo^°\ At the Ri and R2 points, 
there is the contribution to F. At M (0(k) = 18°) the con- 
tribution to phonon softening becomes a maximum. 

Our numerical calculation shows that the RBM of a 
(14,2) chiral SWNT exhibits an E-p dependent frequency 
shift of about 3cm _1 within an E-p range of ±0.3eV. 
This result is consistent with the experimental result by 
Nguyen et alM. who observed that the RBM of a (14,2) 
chiral tube exhibits a small gate dependent frequency up- 
shift up to 3cm _1 within the gate voltage range of ±leV, 
if we assume that the change of the gate voltage of ±1V 
corresponds to E-p ~ ±0.3eV, namely, the gate efficiency 



factor is about 0.3. 



A. Raman Intensity 

We note that Eq. pTj) does not mean that the Raman 
intensity of the RBM is much smaller than that of the 
LO/TO phonon modes. The Raman intensity is relevant 
to the el-ph matrix element for a photo-excited electron 
within the conduction states, 



(ee(k)\n cp \J 0) ) 
= J (^ c K k (r))tW cp V, c K k (r)d 2 r 



2r \ -er 



«-»ik) 



goff „-H38 



goff r -i38 
2 e 




(2g on + g oS cos(6(k) + 30)) . 



<->(ki 



'-'Ik) 



(25) 



In this case, the g D n term does contribute to the matrix 
element and enhances the Raman intensity. Although 
Eq. (|25p is chirality dependent, the dependence is small 
since 2g on ^> g a g. We compare this result with the cor- 
responding LO mode el-ph matrix element from a con- 
duction state to a conduction state in zigzag or armchair 
SWNTs near the K point ^ 



(ee(k)|H op |w LO ) 



The ratio between Eqs. 



\(ee(k)\n c M 0) )\ ffona 



ig s— — cosO(k). (26) 

Ore. 



and 



becomes 



|(ee(k)|W e p|wT J o)| 9oS r u L o 



(27) 



which means that the intensity of the RBM can be com- 
parable to that of the G band because the intensity ratio 
is given by 



^RBM 



(ee(k)|We 



,(0)M 2 



3.5(nm) 



(ee(k)|W cp |w L o) 



(28) 



When 7r-electrons satisfying 8(ku) w contribute to 
the intensity most effectively, the intensity of the RBM is 
maximum for zigzag nanotubes (9 = 0) and is minimum 
for armchair nanotubes (9 = 7r/6). The same tendency 
is observed by the experiment of Strano et ali^ and is 
confirmed by a first-principles calculation by Machon et 
al£^ On the other hand, the frequency shift of the RBM 
is absent for armchair nanotubes while that for zigzag 
nanotubes is a maximum. 



V. EXPERIMENT 

We next compare our theoretical results with experi- 
mental Raman RBM data taken for an individual SWNT 
as a function of the gate voltage with the laser excitation 
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energy 2.15eV. In Fig. [Ha), we show a plot of the RBM 
Raman frequency as a function of the gate voltage. It 
should be noted that the En values seem to be modified 
slightly as a function of the gate voltage, which changes 
the Raman intensity. The experimental details relevant 
to Fig. [6] will be reported elsewhere. 43 



(a) 



246 r 
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(2 
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(b) 
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FIG. 6: (Color online) (a) The RBM frequency and (b) the 
spectral width as a function of the gate voltage for an isolated 
SWNT. The laser excitation energy is 2.15eV. The (n, m) 
value of the SWNT is assigned as (11, 2) or (12, 0). 

Here we note that the Raman frequency in Fig. [Ha) 
shows a minimum value at zero gate voltage, while the 
spectral width in Fig.^b) appears to be constant within 
the error bars. The gate voltage independent offset of the 
spectra width does not originate from the RBM phonon 
self-energy due to the electron-hole pair creation but 
rather from the life-time of a photo-excited carrier. The 
absence of a gate voltage dependent broadening shows 
that £ga P is larger than 7ki/°). This is consistent with the 
fact that the (n, m) for this SWNT is assigned as (11, 2) 
or (12,0) using a Kataura plot based on the extended 
tight-binding models In fact, as shown in Fig. [DJb), the 
calculated E gap (62meV for (11,2), 68meV for (12,0)) 
is larger than the RBM energy (ftu/°) ~ 30meV). The 
RBM frequency difference between the gate voltage of 



0V and of IV is about 2.5 cm 1 , which is close to our 
theoretical estimation for a (11,2) or (12,0) SWNT if 
the gate coupling efficiency factor of this SWNT is about 
0.2. We estimated the gate coupling efficiency to be 
about 0.2 for this sample using Eq. (29) in Ref. I22I 
which is an estimation based the voltage window in 
which the broadens^ 3 - In this case, a frequency shift 
is estimated to be about 2.4cm -1 for a (12,0) SWNT 
from (^ F /0.6eV)50(a cc /d t )cm- 1 with E F = 0.2eV and 
d t ~ 0.93nm. The sample does not exhibit a perfect 
symmetric shape for the RBM frequency about positive 
and negative gate voltage values but shows some asym- 
metric shape. This asymmetry can not be explained by 
the phonon self-energy. We think that the asymmetry is 
due to a change of a spring force constant by doping. 



VI. DISCUSSION AND SUMMARY 

Let us discuss the similarity between the RBM and the 
LO/TO phonon modes for achiral SWNTs. The RBM 
matrix element of Eq. (|23p shows that the 8(k) depen- 
dence is the same as the LO (TO) matrix element of 
Eq. when 9 = 0° (6 = 30°). As a result, we obtain 
using Eq. (JSTJ that 



,(2) 

-'rbm 



(2) O cc 



J LO 



(29) 



for zigzag SWNTs. This correspondence originates from 
the character of the RBM as an optical phonon mode 
like the LO/TO modes; that is, the A-atom and B-atom 
oscillate in the opposite circumferential directions. For 
chiral SWNTs (6 ^ 0° or 30°), on the other hand, the 
relationship between the RBM and the LO /TO is not as 
straightforward as that for achiral SWNTs because the 
LO/TO phonon eigenvector is not pointing along either 
the nanotube axis or the circumference 4^ This changes 
the O(k) dependence of the LO (TO) matrix element of 
Eq. Q3JJ as 



(eh(k)|H C p|w LO ) = -ig D ff sin (6(k) - 

^CC 

(eh(k)|H C p|w T o) = -«3off — ^ cos (9(k) 

OjCC 



(30) 



where is a parameter for the phonon eigenvector^ If 
we consider the frequency shift for \E-p\ ^ 7io;lo(to)/2, 
only high-energy electron-hole pairs contribute to the fre- 
quency shift. Then, by rewriting Eq. (|23[) using Eq. (I30[) . 
and putting the result into Eq. we obtain using 
Eq. that 



,(2) 
RBM 



(2) °cc 2 /on 



(31) 



(21 

where we assume <fi < 30° which results in |wrpo| <C |w£ 
We note that the relationship between the frequency shift 
of the RBM and that of an optical phonon mode, which 
is similar to Eq. (|3"Tj) , was also derived by Nisoli el al^ 
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We have shown that the electron-hole pair creation is 
given by the off-site deformation potential, whose effect 
is represented by the deformation-induced gauge field, 
A(r). It is naturally expected that electron- hole pair cre- 
ation is enhanced where the deformation-induced gauge 
field appears. A static A(r) field appears near the bound- 
ary (edge) of the sampled We will study the effect of the 
edge on the phonon frequency shift in the future. 

Finally we discuss the effect of impurities on the RBM 
frequency. When impurities are approximated by adding 
V(r)ao to the effective-mass Hamiltonian of Eq. (|4|), then 
it can be shown that the broadening of the RBM due to 
the el-ph interaction is not enhanced by the presence of 
impurities. It is because of this that the wavefunction of 
a conduction state relates to the wavefunction of a va- 
lence state by multiplying a z (particle-hole symmetry), 
and therefore the matrix element of V(r)<ro between a 
conduction state and a valence state vanishes. In con- 
trast, the Coulomb interaction between an electron and 
a hole may contribute to the broadening of the RBM 
since an electron and a hole in the intermediate state are 
attracted to each other, which affects the life time of the 
RBM. 

In summary, for a fixed diameter metallic tube, a 
zigzag SWNT exhibits the maximum RBM frequency 
shift and an armchair SWNT does not show any fre- 
quency shift. This is due to the chirality dependent el-ph 
interaction for the RBM. For a zigzag SWNT, the fre- 
quency softening is about 10 cm -1 in a (9,0) SWNT 
and it is proportional to the inverse of the diameter. 
When dt ;$ 2nm, no broadening of the RBM spectra 
appears, since the curvature-induced energy gap is larger 
than the RBM phonon energy. (9, 6) and (10, 7) SWNTs 
are candidates which can exhibit a broadening since the 
curvature-induced energy gap is smaller than the RBM 
energy. Although the frequency shift of the RBM is much 
smaller than that of the T point optical phonon mode, 
the frequency shift of the RBM as a function of the Fermi 
energy shows a characteristic behavior depending on the 
relative position of the cutting line with respect to the 
Dirac point. 



formation, the Hamiltonian of 7r-electrons is modeled by 



Ho 



-7o 2^ 

ieA a=l,2,3 



E (( C ?+«M + (c^cf +a ) . (Al) 



We use the Bloch theorem to diagonalize Eq. (|A1|) . The 
Bloch wavefunction with wavevector k is defined by 




(cf)t|Q) ( S = A,B), (A2) 



where N u is the number of hexagonal unit cells and |0) 
denotes the state of carbon atoms without 7r-electrons. 
The off-site matrix element of 7io is given by 



(n\n \n) = 



(*bIWa> 



where / a (k) 



-7o £ /a(k) = -7o/(k), 

a=l,2,3 

-7o £ 

a=l,2,3 



_ / a (k)* = -7o/(k)*, 

a=l,2,3 

ik R » and /(k) = £ 3 / a (k). 



(A3) 
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and 



the on-site matrix element of H , (^\Ho\^) (s =A,B), 
can be taken as zero. The energy eigenequation is written 
in 2 x 2 matrix form as 



E(k) 



*b! 



-Oo 



/(k)\ 
/(k)* J I |*g 



(A4) 



The energy band structure of H.q is obtained by solving 

0. (A5) 



det ( #(k) 7o/(k)\ 
d6t { 7 of(k)* E(k) ) 

The solution, E{k) = + 7o |/(k)| (- 7o |/(k)|), is the con- 
duction (valence) energy band. The conduction energy 
band and the valence energy band touch each other where 
|/(k)| vanishes. |/(k)| = is satisfied at the K point, k F 
(= (47r/3V3a cc ,0)), and at the K' point, -k F . The K or 
K' point is referred to as the Dirac point. 

By expanding f a (k) in Eq. (|A3|) around the wavevector 
of kp (the K point), we obtain / a (k F + k) = / a (k F ) + 

i/ a (k F )k ■ R a H . Using k F = (47r/3V3a cc , 0), we 

get A(k F ) = 1, / 2 (k F ) = e~ 1 ^, and / 3 (k F ) = . 
Substituting these into Eq. (|A3[) . we obtain 
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APPENDIX A: EFFECTIVE MASS THEORY 
WITH LATTICE DEFORMATION 

In this Appendix we derive Eqs. ([5]) and © from 
the nearest-neighbor tight-binding Hamiltonian of a 
graphene sheet with a lattice deformation. 

When a graphene sheet does not have any lattice de- 



A k |^0|^ F+k ) 



3o<cc / , 
70— 7T- \ k x 



^k F +k 



2 

7o— ^(k x + iky 



(A6) 



0. We ne- 
U) which is 



where we used (4-^ l^ol*^) = "7o/(k F 
gleet the correction indicated by • • • in Eq 
of order of 0(k 2 ) because we only consider k-states near 
the K point, namely |k| <C |k F |. 
From Eq. (|A6[) . we see that Eq, 

by 



H is approximated 



E(k F + k) 

3 7 dec 



^kp+k^ 



k F +k> 
B 



k^£ I z k'u 



(A7) 
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By introducing the Fermi velocity as u F = 3 7 oa cc /2?i, the 
momentum operator p = —ihV, and the Pauli matrix 
v = (^xiO'y), we obtain the effective- mass Hamiltonian 
as upff • p which is given by Eq. 

A lattice deformation induces a local modification of 
the nearest-neighbor hopping integral as —70 — * —70 + 
*7§(ri) (a = 1,2,3) (see Fig.^a)). We define this per- 
turbation as Eq. ([7]). The off-site matrix element of Hi 
with respect to the Bloch wavefunction of Eq. (|A2[) is 
given by 



(^ +4k |Wi|*S) 

= E <57o a (ri)/ a (k)e 

u igA o=l,2,3 

(*5 +4k |«l|*A> 

57o a (ri)/o(k)*e-« k -(-+ R *). 



— i5k-ri 



(A8) 



igA a=l,2,3 



By changing k in Eq. (|A8|) to k F + k and using / a (k F 
k) = / (kp) + i/ a (k F )k • R a H , we see that 



4EE <57 a (rO/ a (k F )e-' ;5k ^ 

u igAo=l,2,3 
^k P +k+5k| Wi |^k F +k^ 

-J-H *7-(ri)/-(kp)*e- iak -" 



(A9) 



igA o=l,2,3 



The correction indicated by • • • in Eq. (IA9|) is negligible 
when |k| < |k F | and |<5k| < |k F |. Substituting A(k F ) = 
1, / 2 (k F ) = e- l2 f and / 3 (k F ) = e +iS r into Eq.JSH]), we 
get 



( ^k F+ k+ 5 k| Wi |^K, 



k F +k\ 



<5kr 



igA 



^k P +k+ik| Wi |^k P +kj 



(A10) 



— i<5k-r 



igA 



where A(r) = (^4 x (r), A y (r)) is defined from Sj§ (r) (a 
1,2,3) as 



VF A x (r) = ^(r) - - (5jZ(v) + £ 7o 3 (r)) , 
v F A y (r) = ^-(6 1 2 (r)-6 1 3 (r)). 



(All) 



Equation (|A10jl shows that Hi appears as vpcr ■ A(r) in 
the effective-mass Hamiltonian. Thus, the total Hamilto- 
nian becomes vpcr ■ (p + A(r)), which is given by Eq. ([5]). 



APPENDIX B: CHIRALITY DEPENDENCE OF 
THE ELECTRON-PHONON INTERACTION 

Next we derive the el-ph interaction for the RBM in 
an (n,m) SWNT. 

A straightforward method to derive the effective-mass 
Hamiltonian for an (ra, m) nanotube is to represent a vec- 
tor in terms of unit vectors ei = C/j/|C^| and e 2 = 
T/|T| instead of and e y (see Fig. EJa)). Here C/j and 
T are the chiral and translational vectors, respectively.— 
The relationship between (ei, e 2 ) and (e^, e y ) is given by 
the chiral angle as 



cost/ — sm( 
sin 9 cos 9 




(b) 



(Bl) 




FIG. 7: (a) We define the coordinate system (£1,0:2) for 
graphene. ei (e2) is the dimensionless unit vector for the x\- 
axis (o;2-axis). (b) x\ (X2) is the coordinate around (along) a 
tube axis. The periodic boundary condition along the x\ axis 
corresponds to a carbon nanotube with chiral angle 6. 

Using Eq. (|B1[) . we represent R a (a = 1, 2, 3) in terms 
of ei, e 2 , and 9 as 



^ • a 1 , 
— — sin — - cos ( 

2 2 

^ • n 1 , 

sin cos ( 

2 2 




(B2) 



Then, by following the same procedure to get Eq. (|A6 
of Appendix [XJ we obtain 



(*j? +k |«o|*l P+k ) = e +ie 7 o^(fci ^- ,/. >) + . 



(B3) 



where we denote the wave vector in the direction around 
(along) the tube axis fcj (fc 2 ) that is k = fcxei + & 2 e 2 (see 
Fig. (Ha) and (b)). 
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Let us introduce a new Bloch wave function which is 
defined by adding a chiral angle dependent phase to the 
original Bloch wave function of Eq. (|A2j) as 



l*A(0)>=e- <fl/ W>, 
I*b(0)> 



= e +ie ' 2 \^l). 
Then the effective-mass Hamiltonian becomes 

v?(a x pi + a y p 2 ), 
and Eq. (|A10[) becomes 

($kF+k+<5k^| Wi |^k F +k^^ 



ieA 



k P +k+5k^| Wi |^k P +k^ 

W, 



ieA 



(B4) 



(B5) 



(B6) 



Therefore, by introducing A\(r) and -42 (r) which are de- 
fined by 



Ai(r)\ _ ( cos 9 sin6>\ ( A x {r) 



Mr) 



sin cos ( 



A y (r) ' 



(B7) 



we see that the off-site interaction can be written as 



Hf = a x A l (r)+a y A 2 (r) 



(B8) 



The Hamiltonian for a chiral SWNT can be represented 
by 



v F a- ■ (p + A(r)) 



where p = (pi,p 2 ) and A(r) = (A^r), A 2 (r)). 



(B9) 



To determine A\ (r) and A 2 (r) , it is necessary to repre- 
sent A x (r) and A y (r) in the new coordinate system. We 
shall do this for obtaining the el-ph interaction for the 
RBM. Putting Eq. (JB2J) into 



^7o a (r)«ffRa-((Ra-V)s(r)) (BIO) 



of Eq. {TTJ , we get using Eq. © 



v F A x (r) 
v F A y (r) 



ffoff 

2 

9oS 



cos 26* 



sin 26* 



dsi(r) ds 2 (r) 
9si(r) 



9a; i 



dx 2 
ds 2 (r 
dx 2 



sin 26» 



dsi(r) <9s 2 (r) 



cos 26» 



dx 2 
gfi(r) 



dxi 
ds 2 (r 
dxi 



(Bll) 



Inserting Eq. (|B11|) into Eq. (|B7j) . we see that 



(_9si(r) ds 2 (r) \ 
fai(r) ds 2 (r) ' ^ j 
9x 2 9xi ' 



The off-site el-ph interaction for the RBM of a chiral 
SWNT is given by replacing 

dsi(r) _^ 0fli(r) + 
dx\ dx\ r 



in Eq. (|B12|) as we have explained in Sec. IIIII For the 
T point RBM, since ds±(r)/dxi = and ds 2 (r)/dxi = 
(i = 1,2), we obtain v F Ai = —g g(s z /d t ) cos 39 and 
v F A 2 — g e(s z /dt) sin 39. This gives the el-ph interaction 
written in Eq. (|22"]) . 
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